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We present numerical support for the hypothesis that macroscopic observables of dense granular 
media can be evaluated from averages over typical blocked configurations: we construct the corre- 
sponding measure for a class of finite-dimensional systems and compare its predictions for various 
observables with the outcome of the out of equilibrium dynamics at large times. We discuss in detail 
the connection with the effective temperatures that appear in out of equilibrium glass theories, as 
well as the relation between our computation and those based on 'inherent str ucture' arguments. A 
short version of this work has appeared in Phys. Rev. Lett. 85, 5034 (2000), [sond-mat/000614C . 
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I. INTRODUCTION 

Granular systems involve many particles, so there 
is a strong motivation to treat them with thermodynamic 
methods. This approach is justified when one is able to 
identify a distribution that is left invariant by the dynam- 
ics (e.g. the microcanonical ensemble), and then assume 
that this distribution will be reached by the system, un- 
der suitable conditions of 'ergodicity'. Unfortunately, be- 
cause energy is lost through internal friction, and gained 
by a non-thermal source such as tapping or shearing, the 
dynamical equations do not leave the microcanonical or 
any other known ensemble invariant. Moreover, just as in 
the case of aging glasses, the compaction dynamics does 
not approach any stationary state on experimental time 
scales. 

Consider a compaction experiment, in which we sub- 
ject a granular system to gentle, periodic tapping. To 
keep the discussion simple, we can assume that there is 
no gravity, and that there is a piston applying a constant 
pressure on the surface. The system compactifies very 
slowly ||, in practice never reaching the most dense, op- 
timal packing. At a given long time, when the system 
has density p(t), we may wish to measure for example 
the fraction of grains that are at relative distance r: the 
structure factor. This quantity being averaged over all 
particles, one can expect it to be a reproducible observ- 
able. However, there is in principle no method to calcu- 
late the structure factor other than solving the dynamics. 

Some years ago Edwards proposed that one could 
reproduce the observables attained dynamically by cal- 
culating the value they take in the usual equilibrium dis- 
tribution at the corresponding volume, energy, etc. but 
restricting the sum to the 'blocked' configurations defined 
as those in which every grain is unable to move. In the 
case of the previous paragraph, we would compute the 
structure factor in all the possible blocked configurations 



of density p(t), and calculate the average. Thus, the only 
input from dynamics would be p(t), apart from which the 
calculation is based on a statistical ensemble. 

This 'Edwards ensemble' leads immediately to the def- 
inition of an entropy (in the glass literature a 'com- 
plexity') S e dwi given by the logarithm of the number 
of blocked configurations of given volume, energy, etc., 
and its corresponding density s e dw = S e d w /N. As- 
sociated with this entropy are state variables such as 
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Recent developments in glass theory, especially those 
related to their out of equilibrium dynamics, have come 
to clarify and support such a hypothesis — at least within 
mean- field models (see below). The present paper ad- 
dresses the natural question of whether Edwards' mea- 
sure gives good results for the compaction of finite dimen- 
sional, non mean-field models. The result is that there 
is a class of models for which this is the case. A shorter 
version of this work has been published in [Q . 

The article is organized as follows: we first discuss the 
nature of the assumptions (Subsection I-A), and the ev- 
idence in support coming from mean-field models (Sub- 
section I-B). In subsection I-C we discuss in detail the 
checks already made with Lennard-Jones glasses, in the 
context of the so-called 'inherent structures', and their 
relation with the present approach. 

In sections II and III we treat two finite-dimensional 
models which reproduce many of the features of glasses 
and granular media, namely the Kob- Andersen (KA) j^] 
and Tctris 10 models. We devise a method to count and 
calculate averages over the blocked configurations, explic- 
itly constructing in this way Edwards' measure. We com- 
pare expectation values thus obtained with equilibrium 
values and with the outcome of slow, aging dynamics and 
find very good agreement between the predictions of Ed- 
wards' measure and aging dynamics. 
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In order to show that this agreement does not hold for 
all forms of slow dynamics, we repeat the procedure in 
section IV for another model exhibiting slow, logarithmic 
relaxations, the random field Ising model (RFIM) JlT| . 
We conclude with a discussion of our results in section 
V. 

A. The Assumption 

One possibility of making an assumption a la Edwards 
would be to consider a fast quench, and then propose that 
the configuration reached has the macroscopic properties 
of the typical blocked configurations. This would imply 
that the system stops at a density for which the number 
of blocked configurations is maximal. We do not follow 
this path, as we will give sufficient evidence that generi- 
cally the vast majority of the blocked configurations are 
much less compact than the one reached dynamically, 
even after abrupt quenches. 

Our strategy here is instead to quench the system to 
a situation of very weak but non-zero tapping, shearing 
or thermal agitation. In this way, the system keeps com- 
pactifying, albeit at a very slow rate. In this context, 
we consider a flat measure over blocked configurations 
conditioned to having the energy and/or density of the 
dynamical situation we wish to reproduce. This means 
that we have given up trying to predict the dynamical 
energy or density by methods other than the dynamics 
itself. 

That configurations with low mobility should be rele- 
vant in a jammed situation is rather evident, the strong 
hypothesis here is that the configurations reached dy- 
namically are the typical ones of given energy and den- 
sity. Had we restricted averages to blocked configura- 
tions having all macroscopic observables coinciding with 
the dynamical ones, the construction would exactly, and 
trivially, reproduce the dynamic results. The fact that 
conditioning averages to the observed energy and den- 
sity suffices to give, at least as an approximation, other 
dynamical observables is highly non-trivial. 

At this point it is important to warn the reader about a 
misconception. It goes like this: The system has at every 
energy and density many blocked configurations. Now, 
we know that in systems with many minima (for example 
from mean-field glasses) all the minima of given energy 
tend to have the same values for macroscopic observables. 
Hence, it is natural to assume that their basins of attrac- 
tion are themselves also the same, and hence Edwards' 
'flat average ' assumption is justified. 

To see the danger of such a reasoning, let us para- 
phrase it in another context, in which it is clear that the 
conclusion is generically erroneous: consider a driven, 
macroscopic system Jl2[ (e.g. fully developed turbulence) 
thermostated at energy E. By the same token, we would 
say: The system is restricted to move in the energy = E 



shell. Now, we know that almost all points in an energy 
shell of a macroscopic system have in the thermodynamic 
limit the same values of macroscopic observables ( we dis- 
regard symmetry breaking cases). Hence, it is natural to 
suppose that the dynamic stationary measure is also the 
same in all points: therefore the stationary distribution 
is flat, i. e. microcanonical. 

This is of course wrong: we know that generically the 
stationary measure of a driven, thermostated system is 
dominated by an ensemble of zero volume within the en- 
ergy shell. Almost all points in the energy shell do have 
the same weight, but this weight is zero. In order that 
the stationary measure coincides with the microcanoni- 
cal measure we need some strongly specific properties for 
the dynamics: such is the case of chaotic Hamiltonian 
dynamics. The same can be said about Edwards' mea- 
sure: that all blocked configurations of a given energy 
have the same basin of attraction may be quite generally 
true, but in order for Edwards' measure to be relevant 
the combined basin of attraction of the typical configu- 
rations should not vanish! 

Starting from a random configuration the probability 
of falling into a basin is proportional to its volume, so 
with a quench we are not sampling a typical basin: rare 
basins (exponentially smaller in quantity) of large (expo- 
nential) volume may dominate — and this is generally 
true for a quench from equilibrium at any temperature. 
Using Edwards' measure is justified when for some reason 
one can consider that typical basins of a given level are 
also typically accessed: a very strong assumption that is 
generally not valid. The reason this hypothesis (if true) is 
useful is that one can construct averages over configura- 
tions defined by a local property (being blocked) without 
having to know the basin of attraction, which involves 
solving the dynamics. 

There is still however a puzzling question: We are using 
blocked configurations as a distribution for the dynamic 
situation. However, this seems odd, since we know that 
neither a relaxing nor a gently driven system will stay 
in one of those configurations! Here the example of the 
stationary measure in dynamic systems is also instruc- 
tive: we know that such a measure can be constructed by 
considering only the periodic trajectories, although the 
probability that the system is in a periodic trajectory is 
strictly zero. Somehow, these trajectories form a 'skele- 
ton' of the true distribution — and such would also be the 
role of the blocked configurations in Edwards' measure. 

Moreover, we will check that configurations with a 
small, though non-zero fraction of mobile particles yield 
the same statistics. 

B. Solvable models 

As mentioned above, the fact that dynamically ac- 
cessed blocked configurations are the typical ones does 
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not follow from any general principle that we know, and, 
as we shall see below, is indeed not always true. 

In order to progress, one can exploit the analogy be- 
tween the settling of grains and powders, and the aging 
of glassy systems [0 since in both cases, the system re- 
mains out of equilibrium on all accessible time-scales, and 
displays very slow relaxations. 

In the late eighties, Kirkpatrick et al. [ p"3| , |i"4| recognized 
that a class of mean-field models contains, although in a 
rather schematic way, the essentials of glassy phenomena. 
When the aging dynamics of these systems was solved 
analytically a feature that emerged was the existence of 
a temperature Td yn for all the slow modes (corresponding 
to structural rearrangements) fl5|,|l6| . 

For the purposes of this paper, Td yn can be defined 
by comparing the random diffusion and the mobility 
between two widely separated times t and t w of any 
particle or tracer in the aging glass. Surprisingly one 
finds in all cases an Einstein relation ((r(t) — r(t w )) 2 ) = 

Td y n S ^Jf^ tw ^ ; where r is the position of the particle 
and / is a constant perturbing field, and the brackets 
denote average over realizations. While in an equilib- 
rium system the fluctuation-dissipation theorem guaran- 
tees that the role of Td yn is played by the thermodynamic 
temperature, the appearance of such a quantity out of 
equilibrium is by no means obvious. Td yn is different 
from the external temperature, but it can be shown to 
have all other properties defining a true temperature [fl6| . 

As it turned out, despite its very different origin, this 
temperature matches exactly Edwards' ideas. One can 
identify in mean-field models all the energy minima (the 
blocked configurations in a gradient descent dynamics), 
and calculate l/T e d w as the derivative of the logarithm 
of their number with respect to energy. An explicit com- 
putation shows that T e d w coincides with Td yn obtained 
from the out of equilibrium dynamics of the models aging 
in contact with an almost zero temperature bath [jl?] . 
Moreover, given the energy E{t) at long times, the value 
of any other macroscopic observable is also given by 
the flat average over all blocked configurations of energy 
E(t). Within the same approximation, one can also treat 
systems that like granular matter present a non-linear 
friction and different kinds of energy input, and the con- 
clusions remain the same |^3| despite the fact that there 
is no thermal bath temperature. 

Edwards' scenario then happens to be correct within 
mean-field schemes and for very weak vibration or forc- 
ing. The problem that remains is to what extent it car- 
ries through to more realistic models. In this direction, 
there have been recently studies 0] of Lennard- Jones 
glass formers from the perspective of the so-called 'in- 
herent structures' (2j|, which suggest that whatever the 
measure for the slow dynamics, it is not sensitive to the 
details of the thermal history -see next subsection. 

The path we follow here jl]] is to construct the Ed- 



wards measure explicitly in the case of representative 
(non mean-field) systems, together with the correspond- 
ing entropy and expectation values of observables. We 
thus obtain results that are clearly different from the 
equilibrium ones, and we can compare both sets with 
those of the irreversible compaction dynamics. 

C. Differences and similarities with approaches 
based on 'inherent structures' 

Let us discuss in detail the relation between the present 
approach H with the one followed by Kob et al. E3] 
with Lennard-Jones glasses, later applied in the granular 
matter context in |26| . 

One can describe the work of Kob et al. as follows: 
starting from an equilibrated system at temperature T 
above the glass transition, the system is first quenched 
to a temperature 7/ below the glass transition, where the 
system spends aging a time t w , after which it is quenched 
again to zero temperature (In some of the procedures 
there is no intermediate stop: t w = 0). The dynamics of 
the final quench being at zero temperature, the system 
eventually lands in a blocked configuration of energy E. 
For each cooling protocol parameterized by (T,Tt,t w ) 
the final energy £^,27, t„ is a reproducible quantity in 
the thermodynamic limit. 

Suppose now we classify all thermal histories accord- 
ing to the energy E of the blocked configuration reached 
at the end. Kob et al. then ask the following question: 
are all other macroscopic observables fully determined 
by E, or, otherwise stated, is the effect of the whole 
history (T,Tf,t w ) completely encoded in El For the 
macroscopic observables they considered, their answer is 
within numerical precision affirmative. (In their work, 
they chose as macroscopic observables the spectrum of 
the energy Hessian, instead of the structure factor as we 
do here). 

We can now discuss the relation of their approach and 
the present paper. On the one hand, because there is 
no direct sampling of typical configurations of energy 
E, but a comparison of configurations reached after dif- 
ferent histories, the procedure of Ref. does not ad- 
dress the question of what the actual distribution is for 
given E (unless, of course, one makes extra assumptions). 
One could imagine a situation in which a small subset of 
blocked configurations of given energy contributes to the 
measure because they have a larger basin of attraction 
for every history. Indeed, we shall see later that the dy- 
namics of the RFIM, while passing the test of Kob et 
al., is not well reproduced by a flat Edwards' measure. 
What Ref. |24|] does suggest is that whatever the mea- 
sure, it is insensitive to the details of the thermal history, 
which only has the effect of specifying E — two thermal 
histories finishing in the same energy E yield the same 
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values for all the other macroscopic observables. The ap- 
proaches are clearly complementary: the suggestion of 
the present work that Edwards' measure gives good re- 
sults for a slow compaction would be of little use without 
the insensitivity of the measure on the history suggested 
in Rcf. @. 

Let us add that Kob et al. define a temperature T* 
for a process (T,Tf,t w ) by demanding that E obtained 
by a direct quench from T* to zero temperature be equal 
to £r,Tf,t TO - This temperature is not equal to (but may 
be an approximation of) the Edwards temperature which 
we calculate below. 



II. KA MODEL 

The first model we consider is the so-called Kob- 
Andersen (KA) model || that was first studied in the 
context of Mode-Coupling theories |2lj as a finite dimen- 
sional model exhibiting a divergence of the relaxation 
time at a finite value of the control parameter (here the 
density); this divergence is due to the presence in this 
model of the formation of "cages" around particles at 
high density (the model was indeed devised to reproduce 
the cage effect existing in super-cooled liquids). 

Though very schematic, it has then been shown to re- 
produce rather well several aspects of glasses Q, like 
the aging behaviour with violation of FDT j25|, and of 
granular compaction [ |30| . 

The simplicity of its definition and the fact that it is 
non mean-field makes it a very good candidate to test Ed- 
wards' ideas: in fact, the triviality of its Gibbs measure 
will allow us to compare the numerical data obtained for 
the dynamics and for Edwards' measure with the analytic 
results for equilibrium. 



A. Definition 

The model is defined as a lattice gas on a three dimen- 
sional lattice, with at most one particle per site. The 
dynamical rule is as follows: a particle can move to a 
neighboring empty site, only if it has strictly less than m 
neighbours in the initial and in the final position. Fol- 
lowing H, we take m = 4: this ensures that the system 
is still ergodic at low densities, while displaying a sharp 
increase in relaxation times at a density well below 1. 
The dynamic rule guarantees that the equilibrium distri- 
bution is trivially simple since all the configurations of 
a given density are equally probable: the Hamiltonian is 
just since no static interaction exists. 

In order to mimic a compaction (or aging) process 
without gravity, we simulate a 'piston' by creating and 
destroying particles only on the topmost layer (of a cubic 
lattice of linear size L) with a chemical potential p [ p8[ . 
More precisely, each Monte-Carlo sweep is divided in the 



following steps: (i) for each of the site of the topmost 
layer, add a particle if the site is empty, and, if it is occu- 
pied, withdraw the particle with probability exp(— /3/z). 
(ii) try to move each particle, in random order, according 
to the dynamical rule. 



B. Gibbs measure 

Since the Hamiltonian is 0, the equilibrium (or Gibbs) 
measure corresponds simply to a fiat measure over all 
configurations, without taking into account the dynam- 
ical constraint. Therefore, the relation between density 
and chemical potential is 



1 



1 + exp(-/3/i) 



(1) 



and the exact equilibrium entropy density per particle 
reads 



Sequii(p) = -plnp - (1 - p)ln(l - p) 



with in particular 



equil 



(If) 



= -/3/x. 



(2) 



(3) 



In this model, the temperature 1/(3 is irrelevant since it 
appears only as a factor of the chemical potential and we 
can set it to one throughout. 

Besides, the equilibrium structure factor defined as the 
probability that two sites at distance r are both occupied 
is easily seen to be a constant 



9equii (r) = p 2 



r > 



(4) 



No correlations appear since the configurations are gener- 
ated by putting particles at random on the lattice. It will 
therefore be easy, as already mentioned, to compare small 
deviations from g e quii{r), a notoriously difficult task to 
do in glassy systems. Note that it is also easy to numer- 
ically sample Gibbs' measure at any given density, by 
simply generating at random configurations with fixed 
number of particles. 



C. Non-equilibrium dynamics 

The previously described Monte-Carlo procedure al- 
lows to produce equilibrium configurations, even if the 
dynamical constraint is enforced, as long as p is low 
enough. However, at densities close to p g (~ 0.88), the 
particle diffusion becomes extremely slow due to the ki- 
netic constraints. In fact, the diffusion coefficient is well 
approximated by 



D{p) ~ (p s - PY 



(5) 
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with cf) ~ 3.1 The equilibrium with p > p g is not 
reached by compaction after extremely long times: if a 
chemical potential p such that p(p) > p g is applied, the 
system falls out of equilibrium. Moreover, the density 
obtained at long times depends on the history: the slower 
the chemical potential is raised, the denser the system 
becomes |2c| . 

We are here interested in this out of equilibrium dy- 
namics: we therefore perform a compression, starting 
from low density, by raising the chemical potential up 
to a high value p — 3. Since the equilibrium density at 
p = 3 is much larger than the jamming density p g , aging 
and very slow compaction ensue. We record the den- 
sity p(t), the density of mobile particles p m (i), and the 
spatial structure function gd yn (r, t) defined as the proba- 
bility that two sites at distance r are occupied. Since we 
work at finite sizes, and since particles are always added 
at the same layer, density heterogeneities do appear be- 
tween the topmost layer and the rest of the box if the 
compression is too fast. To avoid any systematic error, 
we use a slow compression (At — 10 4 MC sweeps for 
each increase of Ap = 0.01), and we measure the various 
quantities in the center of the box only, where we checked 
that the system is indeed homogeneous. 

We show in Fig. [l] the parametric plot p m (t) versus 
p(t); at short times and low density, it follows the equilib- 
rium curve (obtained by generating at random configura- 
tions of density p and thus measuring the mean density of 
mobile particles for the Gibbs measure); at low p indeed, 
the relaxation time of the system is smaller than the rate 
of increase of p, so the system has time to equilibrate. 
As the density approaches p g however, compaction slows 
down and p m (t) gets smaller. At large times the system 
is approaching (p ~ p g , p m ~ 0). 
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FIG. 1. Parametric plot of the density of mobile particles 
versus the density, for the equilibrium measure (data obtained 
by generating at random configurations of a given density and 
measuring the number of mobile particles), and during two 
compaction procedures (with Ap/At — 10 -5 and 10 -6 ). 
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FIG. 2. Dynamic structure function obtained in a very slow 
compression, at different times, i.e. different values of p; at 
p = 1.5 the system is still at equilibrium and the correlation 
function is therefore equal to g eq uu , but stronger and stronger 
deviations are observed as p is raised. 

We have also measured the dynamical structure func- 
tion gd yn (r,t), displayed in Fig. ||; induced and spon- 
taneous displacements were measured and compared 
in Q , and these data will be displayed and used in sec- 
tion III-F. 



D. The auxiliary model 

We introduce an 'auxiliary model' which will allow 
us to define the Edwards measure for the KA model. 
In this model particles have energy equal to one if the 
dynamic rule of the KA model allows them to move, 
and to zero otherwise. The Hamiltonian is therefore 
highly complicated, involving next-nearest neighbor in- 
teractions. We can however introduce an auxiliary tem- 
perature l/Paux associated to the auxiliary energy E aux 
(equal to the number of particles that are able to move) 
and perform a simulated annealing, at fixed number of 
particles: at low j3 aU x ah configurations are sampled 
uniformly, while, as {3 a ux grows, the sampling is re- 
stricted to configurations with vanishing fraction of mov- 
ing particles. The Monte Carlo procedure uses non-local 
moves (accepted with a standard Metropolis probability 
min {1, exp(— (3 aux AE aux )}) which allow for an efficient 
sampling: these non-local moves have nothing to do with 
the true dynamics of the original model, and therefore 
the auxiliary model is not glassy. 



5 



0.03 'x 



Cl 0.02 



CO. 



<D 0.01 




b) 



CO. 



0.48 



0.44 



0.4 



0.36 



0.32 



FIG. 3. Thermodynamical properties of the auxiliary 
model, a) Energy per particle e aux vs inverse temperature 
Paux at different particle density p. b) Entropy per parti- 
cle Saux as obtained from thermodynamic integration of the 
energy data. 

In this way, we obtain the equilibrium energy density of 
the auxiliary model, e aux (/3 aux , p) and its entropy density 
Saux{Paux-,p) by thermodynamic integration: 



? (Paux 5 P) 



Sequil (p) ~t~ Paux &aux (ftaux 3 P) 

eaux(f3'aux,P)df3' aux ■ (6) 





where we set 



^equil 



(p) 



(7) 



since the limit [3 aux — ► corresponds to the equilibrium 
measure. In fig. ^ we show a subset of data concerning the 
energy and entropy densities of the auxiliary model. (The 
energy has been computed in the range p G [0.65, 0.95] 
with a step in density Ap — 0.005 and for f3 aux S [0, 20] 
with step A(3aux = 0.1). 

We also evaluate the structure function of the auxiliary 
model g a ux(r, f3 a ux) which is shown in fig. ^ for a density 



p = 0.87. It is clear here as well as in fig. ||that the limit 
Paux — * oo considered in the next subsection is already 
approached for f3 aux ~ 5. 
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FIG. 4. Structure function obtained in the auxiliary model 
at different values of the temperature (3 a ux and for a density 
p = 0.87. The data for /3 aU x = 5, 10 and 20 are indistinguish- 
able. 



E. Edwards measure 

To evaluate the observables with Edwards' measure, 
i.e. the set of configurations where all particles are un- 
able to move, we consider now the limit (3 aux — > oo of 
the observables computed in the auxiliary model. For 
example, the Edwards entropy is then obtained as 

Sedw(p) = Urn s((3 aux ,p) 



Sequil (p) 



(P'aux,P)dP 

an 



since 



0a 



lilll ftaux &aux(,ftaux j P) ■ 



(8) 



(9) 



In fig. |]we plot the Edwards and the equilibrium entropy 
as a function of the particle density. 
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FIG. 5. Edwards entropy per particle of the Kob- Andersen 
model vs. density (full curve). For comparison we also show 
the equilibrium entropy (dashed curve). At high enough 
density the curves are indistinguishable, and join exactly 
only at p = 1. The slope of the tangent to s e dw(p) for 
a generic p allows to extract T e d w (p) from the relation 
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dp 



Comparison of Fig. || with Fig. [I] shows that the most 
typical blocked configurations (p ~ 0.75) are irrelevant 
as far as the compaction dynamics is concerned. 



Similarly, the Edwards measure structure function, 
gedwiy), is obtained as 

9edw(r) = lira g aux (r, f3 aux ) . (12) 

and displayed in fig.|| for various densities. 
Two remarks are in order 

• as the density approaches 1, almost all particles 
become blocked even for Gibbs' measure; more 
precisely, lim^i p^ ibbs = 0; thus, Edwards' and 
Gibbs' measures get closer, which is seen in fig. || by 
the fact that the curves for s equ u and s e d w get very 
close, and lim p _,.i T edw = T Gibbs = 1. Similarly, 
gedw(r) deviates less from g equ u as p increases, as 
is shown in Fig. ||. 

• Edwards' measure is precisely defined as a sam- 
pling over configurations with vanishing fraction of 
mobile particles. We mention here for complete- 
ness the straightforward generalization of Edwards' 
measures as the set of configurations with fraction 
of mobile particles smaller than e (cfr. the quasi- 
states of pfjfl); we then can use the knowledge of 
e-aux {Paux > p) to define (3 e as the value of the aux- 
iliary temperature such that e aux (f3 e ,p) = e, and 
thus define 
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FIG. 6. Edwards structure function g e dw(f) obtained as 
limit of the g a ux(r, (3 aux ) for f3 aU x — > oo. As the density in- 
creases the deviation from g eq uii gets less pronounced. 

Since the relation between chemical potential, temper- 
ature and entropy density at equilibrium is given by (||), 
the natural definition for Edwards' temperature is 



t: 



1 d,S e dw (P) 

p dp 



(10) 



However we work here at fixed density for Edwards' mea- 
sure, and we therefore compute: 



dp 



(ii) 



s ldw(p) = s (Paux =/3 e ,p) 
= SequU (p) +£j3 t 

rP' 

((3aux,p)d(3 

(10 



It is clear that s\ dw 



= > equi l 



a, while s° ed 



(13) 
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We can also measure the structure factors g e (r) 

g(r,P e )- 



F. Comparing the measures 

We are now in a position to compare the long-time 
results of the out of equilibrium dynamics with those 
obtained with the different measures. Section II-C has 
already made clear that the equilibrium measure is not 
able to describe these results. Fig. shows a plot of the 
mobility 



X(t,tw) 



i 



r a k (t w ))) 



a=l fc=l 



Sf 



(14) 



obtained by the application of random forces to the par- 
ticles (see [^| for details), vs. the mean square displace- 
ment 



B{t, t w ) 



1 

3N 



N 



a=l fe=l 
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testing in the compaction data the existence of a dynam- 
ical temperature T c i yn [§9| . (N is the number of particles 
and a runs over the spatial dimensions) One first remarks 
the existence of a dynamical temperature Td yn - Further- 
more the agreement between Td yn and T e d w , obtained 
from the blocked configurations as in Fig. |5[ for the den- 
sity at which the dynamical measurement were made, is 
clearly excellent. 

In Fig. U we plot the long-time dynamical gdyn(r-)t), 
the equilibrium g eq uii(r) — p 2 , and the Edwards' g e d,w{i~) 
structure factors, for the same density p ~ 0.87. While 
gequii(r) is flat, the system has developed during its dy- 
namical evolution some structures, which seem to be re- 
produced rather well by g e dw{r)- Note that using the gen- 
eralized Edwards' measures does not improve the agree- 
ment in a clear-cut way because of the rather large error- 
bars on the dynamical data. 

To summarize, during the compaction, the system falls 
out of equilibrium at high density, and is therefore no 
more described by the equilibrium measure. It turns out 
that Edwards' measure, constructed by a flat sampling 
of the blocked configurations at the dynamically reached 
density, reproduces the physical quantities measured at 
large times, and in particular predicts the correct value 
for the dynamical temperature. 
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FIG. 7. Einstein relation in the Kob- Andersen model: plot 
of the mobility x{t,t w ) vs. the mean-square displacement 
B(t,t-w) (data shown as circles). The slope of the full straight 
line corresponds to the equilibrium temperature (T = 1), and 
the slope of the dashed one to Edwards' prescription obtained 
from figure || at p(t w ) = 0.848. 
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FIG. 8. Structure functions g(r) — p at density p 
computed with the equilibrium, Edwards' and dynamical 
measure of the Kob- Andersen model. The three sets of data 
come from independent Monte-Carlo simulations. The dy- 
namic structure function (circles) is obtained after a very slow 
compression by raising the chemical potential from fi = 1 to 
H = 3 with an annealing rate of 10 -6 Monte Carlo sweeps. 
The Edwards' structure function (open squares) is obtained 
from the auxiliary model. Although the equilibrium value of 
g(r)—p 2 is exactly 0, we also obtain it by a Monte-Carlo simu- 
lation (full squares) in order to show that the difference in the 
short distance behaviour is not an artifact of the numerical 
simulation). The size of the typical error bar on dynamical 
data is shown at r — 3. 



III. TETRIS MODEL 

In this section we extend the results obtained for the 
KA model to another class of models, the so-called Tetris 
Model (TM) p0| . We proceed as before by construct- 
ing explicitly Edwards' and Gibbs' measures, together 
with the corresponding entropy and expectation values 
of some observables, and comparing both sets of data 
with those obtained with an irreversible compaction dy- 
namics. Notice that in this case the equilibrium measure 
is by no means trivial and it has to be computed numeri- 
cally, using an auxiliary model as for the construction of 
Edwards' measure. 



A. Model Definition 

The essential ingredient of the TM |[o| is the geomet- 
rical frustration that for instance in granular packings 
is due to excluded volume effects arising from different 
shapes of the particles. This geometrical feature is cap- 
tured in this class of lattice models where all the basic 
properties are brought by the particles and no assump- 
tions are made on the environment (lattice). The inter- 
actions are not spatially quenched but are determined 
in a self-consistent way by the local arrangements of the 
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particles. It is worth to notice how in this class of mod- 
els the origins of the randomness and of the frustration 
coincide because both are given in terms of the particle 
properties. 

Despite the simplicity of their definition, these systems 
are able to reproduce many general features of granular 
media: the very slow density compaction |1C| ], segrega- 
tion phenomena plfl, dilatancy properties ]33| ] as well as 
memory Q and aging p^JOL 

Let us recall briefly the definition of the model, which 
includes, like in the real computer game Tetris, a rich 
variety of shapes and sizes. On a lattice each particle 
can be schematized in general as a cross with 4 arms (in 
general the number of arms is equals to the coordination 
number of the lattice) of different lengths, chosen in a 
random way. An example of particle configuration on a 
square lattice is shown in Fig. ^. 




FIG. 9. Sketch of a local arrangement of particles in the 
Tetris Model: each particle can be schematized in general as 
a cross with 4 arms of different lengths, chosen in a random 
way. 

The interactions among the particles obey the general 
rule that one cannot have superpositions. For instance 
one has to check that for two nearest-neighbor particles 
the sum of the arms oriented along the bond connecting 
the two particles is smaller than the bond length. It 
turns out that in this way the interactions between the 
particles are not fixed once for all but they depend on 
the complexity of the spatial configuration. 

The extreme generality of the model definition allows 
a large variety of choices for the particles. While the 
original model deals with simple rods, the fact that these 
rods can arrange in an antiferromagnetic-like configura- 
tion of density 1 has motivated the use of random shapes 
to avoid this pathology. In this study, on the other side, 
we will use the so-called "T" -shaped particles defined in 
such a way that three arms have length equal to jd and 
the fourth one zero length, d sets the bond size on the 
square lattice. With this definition one has four types of 
particles corresponding to the four possible orientation 



of the "T"s on a square lattice. Our choice has the fol- 
lowing advantages: on the one hand, no averaging over 
the disorder is needed; on the other hand, the process 
of applying a chemical potential together for particles of 
various sizes could produce a 'filtering' effect that would 
dynamically lead to an artificially dense system with only 
small particles. With the above given rules one can de- 
fine the allowed configurations. One can easily realize 
how the maximal allowed density is equal to p max = | 
which corresponds to a number of possible configurations 
proportional to the linear size of the lattice, L (due to 
the possible translations and symmetries). Fig. |l^ shows 
two possible configurations with p = p rnax . 




FIG. 10. Two possible configurations of the Tetris model 
considered with the maximal density p = p max = 2/3. 

B. Equilibrium measure 

The equilibrium measure is obtained with an anneal- 
ing procedure. We can introduce a temperature T = 1/(3 
associated with an energy E defined as the total parti- 
cle overlaps existing in a certain configuration. For each 
value of T one allows the configurations with a proba- 
bility given by e~ /3E . Starting with a large temperature 
T (a very small (3) one samples the allowed configura- 
tions by progressively decreasing T (increasing f3). As T 
is reduced E decreases and only at T — (no violation 
of constraints allowed) the energy is precisely zero. The 
exploration of the configuration space can be performed 
in two ways: 

1) working at constant density by interchanging the 
positions of couples of particles. This procedure is 
used to compute E(/3,p) and e(/3, p) (energy per 
particle) , from which one can compute the equilib- 
rium entropy per particle by the expression: 

Sequil(p) = S equ u((3 = 00,p) = 

poo 

= s equi i((3 = 0,p)- / e(P,p)d/3 (16) 
Jo 

For the choice made for the particles one has 
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«G8 = 0,P) 



—pin p - 
+pln4, 



(l-p)ln(l-p) 



(17) 



which is easily obtained by counting the number of 
ways in which one can arrange pL 2 particles of four 
different types on N — L 2 sites; 

2) working at constant chemical potential p by adding 
and removing particles. In this case the density 
fluctuates and one measures E((3,n). 

In both cases one can measure the particle-particle cor- 
relation function (equivalent to the void-void correlation 
function) g e quil(f)- Since, working at constant p the den- 
sity fluctuates, we have always measured directly the en- 
semble average < (g e quti( r ) ~ P 2 ) > m order to make 
the data obtained with the two different methods com- 
parable. One can then compare the correlation function 
obtained at constant p (which corresponds to a certain 
average density) with the correlation function obtained 
working at constant density. The results of both methods 
are equivalent. 

Fig.s (|TT| ) and ( |l2| ) report the results for e(/3,p) and 
Sequii(p) as obtained from (|l6|). 
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FIG. 11. e(f3 — 1/T, p) for various values of p. From top to 
bottom the density decreases monotonically from p = 0.6 to 
p = 0.1 
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FIG. 12. Entropy per particle in equilibrium s equ u(p) and, 
for reference, entropy per particle at T aU x = oo. s equi i(p) goes 
to zero at p = |. 



Fig. [13] reports the results for < {g eq uii{i') — p 2 ) > for 
different values of p. The correlation functions g (g e q U u 
or g e dw or gdyn) actually display oscillations around p 2 , 
whose origin can be easily understood: if a particle occu- 
pies a site, the exclusion rules decreases the probability 
that a neighboring site will be occupied. We therefore 
plot in the figures < (g(r) — p 2 ) > 2 in order to show the 
exponential decay of the correlations. 
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FIG. 13. < (g eq uii(r) — p 2 ) > 2 for various values of p. 
From top to bottom the density decreases monotonically from 
p = 0.64 to p = 0.1. 
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C. Edwards' measure 



Edwards' measure is obtained with an annealing pro- 
cedure at fixed density. This means that one samples 
the configurational space by interchanging the positions 
of couples of particles without violations of constraints. 
In this way one is sampling the configurational space 
corresponding already to T = 0. In order to select 
only the subset of configurations contributing to the Ed- 
wards' measure we introduce an auxiliary temperature 
T aU x (and the corresponding (3 aux = l/T aux ) and, as- 
sociated to it, an auxiliary energy E aux which, for each 
configuration, is equal to the number of mobile parti- 
cles, in the same way as for the KA model. A particle 
is defined as mobile if it can be moved according to the 
dynamic rules of the original model. 

Let us describe in detail how the measurements are 
performed. Since for each given density one is interested 
in the subset of the equilibrium configurations with a re- 
duced particle mobility, we start with an annealing pro- 
cedure precisely identical to the one used for the equilib- 
rium measure. This procedure allows us to reach a start- 
ing configuration with a given density and no constraints 
violated. At this stage we perform a Monte Carlo proce- 
dure which exchanges the positions of couples of particles 
without violation of constraints: this procedure accepts 
the non-local moves with a standard Metropolis probabil- 
ity min(l, exp(— f3 a uxAE aux )) which allow for an efficient 
sampling. These non-local moves define an auxiliary dy- 
namics which has nothing to do with the true dynamics 
of the original model, and therefore the auxiliary model 
is not glassy. 

As for the equilibrium we have performed a certain set 
of measures. 

In particular we have measured E aux (f3 aux , p), i.e. the 
decrease of the auxiliary energy at fixed density. In order 
to do this one performs an annealing procedure increasing 
progressively f3 aux and monitoring for each (3 aux the cor- 
responding configurational energy E aux (/3 aux , p) . From 
this measure one can compute the Edwards' entropy per 
particle defined by: 

/•oo 

SequilyP) / &aux(Pauxi P)dj3 aux (1$) 

Jo 

where e aux ((3 aU x, p) is the auxiliary Edwards' energy per 
particle and we put s equ a{p) = s aux {(3 aux =0,p). 

Fig. |TJ reports the results for s e d w (p) as obtained from 
(|l|) compared with s equil (p). 




0.2 0.4 0.6 0.8 

P 

FIG. 14. Edwards' entropy per particle, s e dw{p) and equi- 
librium entropy per particle, s equ u(p). Both go to zero at 

P = l 

For the computation of the particle-particle correla- 
tion function we have to use a different strategy. Always 
starting from a configuration with a given density and 
no constraints violated one performs a Monte Carlo pro- 
cedure, at (3 aux fixed, which exchanges the positions of 
pairs of particles without violation of constraints. Each 
single simulation uses p and f3 aux as input parameters. In 
practice one is trying to sample all the configuration of 
density p and with a particle mobility defined by I3 aux . In 
this context the Edwards' prescription should correspond 
to the limit (5 aux — * oo. In this way we have computed 
< {gedw{r) — p 2 ) > for several values of /3 aux and we re- 
port the results in Fig. [l5]. As one can see the limit 
Paux — > co is attained already for aux of the order of 6. 
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FIG. 15. < {g edw (r)-p 2 ) > 2 for p = 0.58; (3 aux = 1,2,4,6. 
The behaviour seems to saturate to a limit that should cor- 
respond to f3 aux — > oo already for p aux of the order of 6. 

D. Irreversible Compaction Dynamics 

We now turn to the out of equilibrium dynamics of 
compaction, simulated by starting from an empty lattice 
and performing N steps of attempted particle additions 
followed by M steps of attempted particle diffusions. Fig. 
|i~6| reports the results for the density increase and for the 
fraction of mobiles particles as a function of the density 
for N = 1, M = 1 and N = 10, M = 1. While the limit 
N = fixed M — > oo should coincide with the equilib- 
rium, it is clear that the system, at finite M, is not able 
to approach its maximum density and falls out of equi- 
librium. Here again the density at which the number of 
blocked configurations is maximal is smaller than the one 
achieved by compaction. 



FIG. 16. Pmob vs. p for two irreversible dynamics with 
N = 1, M = 1 and N = 10, M = 1. For reference is plotted 
the equilibrium curve. 

It is particularly interesting to notice that in the out-of- 
equilibrium configurations visited during the irreversible 
dynamics the fraction of mobile particles p mo b at fixed 
density is systematically smaller than the corresponding 
value in equilibrium. This suggests the possibility of dis- 
tinguishing between equilibrium and out-of-equilibrium 
configurations by looking at the spatial organization of 
the particles in both cases. We have then measured, dur- 
ing the compaction dynamics, the particle-particle corre- 
lation function at fixed density. In the next section we 
shall compare these results with those obtained on the 
basis of the equilibrium and Edwards' measures. 

E. Comparing different measures 

We are now able, as for the KA model, to compare 
Gibbs' and Edwards' measures with the results of the 
out-of-equilibrium dynamics at large times. 

In Fig. [I?] we plot the deviations of the particle-particle 
correlation functions from the uncorrelated value p 2 . In 
particular we compare < (gdyn{r) — p 2 ) > 2 obtained dur- 
ing the irreversible compaction (N = 1, M = 1) with the 
corresponding functions obtained with the equilibrium 
and Edwards' measures. It is evident that the correla- 
tion function, as measured during the irreversible com- 
paction dynamics, is significantly different from the one 
obtained with the equilibrium measure. On the other 
hand the correlation functions obtained with Edwards' 
measure are able to better describe what happens during 
the irreversible dynamics. In particular what is observed 
is that the correlation length seems to be smaller for con- 
figurations explored by the irreversible dynamics than in 
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the equilibrium configurations. This aspect is captured 
by Edwards' measure which selects, the better the larger 
Paux, configurations with a reduced particle mobility. In 
practice one can summarize the problem as follows: given 
a certain density, one can arrange the particles in differ- 
ent ways. The different configurations obtained in this 
way differ in the particle mobility and this feature is re- 
flected by the change in the particle-particle correlation 
properties. 
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FIG. 17. Comparison between the correlation functions ob- 
tained with the equilibrium measure, the Edwards measure 
(/3aux = 6) and the irreversible dynamics (N = 1,M = 1). In 
all cases the system is considered at a density of p ~ 0.58. 

Also in this case, as in the KA example, it turns out 
that Edwards' measure, constructed by a flat sampling of 
the blocked configurations, is able to reproduce the phys- 
ical quantities measured at large times. Investigations 
are currently running to check also in the Tetris model 
whether the temperature predicted with the Edwards' 
approach coincide with the dynamical temperature Td yn 
as defined for the KA model. 

IV. THE 3-DIMENSIONAL RANDOM FIELD 
ISING MODEL AT LOW TEMPERATURE 

In this section, we consider a case in which Edwards' 
ensemble does not give good results: the low tempera- 
ture domain growth dynamics of a 3D Ising model in a 
weak random magnetic field. This model has been ap- 
plied in many different contexts [ [ill ; and in particular in 
relation with glasses |36) . The logarithmic relaxations it 
displays at low temperature, as well as the dependence on 
its thermal history (like the influence of the cooling rate) 
can also induce comparisons with granular compaction. 



The model is defined as usual: N = L 3 Ising spins 
(si = ±1) sitting on the sites of a regular lattice of linear 
size L, interact ferromagnetically, in a random external 
field. The Hamiltonian is 

H = - J 2J SiSj -'^hiSi , (19) 

<i,j> i 

where the sum runs over pairs of nearest neighbours. 

The strength of the ferromagnetic interaction J can be 
set to 1, and the distribution of the random fields hi will 
be taken as bimodal hi = ±/io- At high temperature, 
the system is in a paramagnetic phase; at low tempera- 
ture and weak magnetic field, there exists a ferromagnetic 
phase (see |ll]]). The typical equilibrium configurations 
are therefore magnetized. 

In the absence of random fields, the low temperature 
dynamics is the well known coarsening of domains of plus 
or minus spins, whose typical size grows as a power of 
time. In the case of the RFIM, the domain walls are 
pinned by the field, and the dynamics proceeds by ther- 
mal activation. The size of the domains therefore grows 
only logarithmically with time. Moreover, the fact that 
thermal barriers are easier to overcome at not too low 
temperatures induces a strong dependence on the cool- 
ing rate. As above, we are interested in the limit of low 
but non-zero temperatures. 

In a large system, the long-time configurations ob- 
tained dynamically are intertwined domains of 'up' and 
'down' spins having similar volumes, the global magneti- 
zation being zero. This is quite different from the equilib- 
rium configurations at the same energy, which are mag- 
netized. In fact, an easy way to show that the long-time 
dynamical configurations differ from the equilibrium ones 
is to look at the distribution of the magnetizations P{M) 
in both cases: at equilibrium P equ u{M) is peaked around 
±M(T), with M(T) > 0, while for the dynamics one ob- 
tains at any finite time for Pd yn (M) a single peak around 
M = 0: in this domain growth dynamics, the system does 
not choose at any finite time between the two basins of 
attractions of the two ground states |}7) . 

The dynamics proceed by thermal activation; therefore 
the long-time configurations are 'blocked', in the sense 
that, at zero temperature, the system would be unable 
to escape from them. The question in the present context 
is now whether these 'blocked configurations' are typical, 
i.e. if their characteristics are well reproduced by a flat 
sampling of all blocked configurations of the same energy. 
We have therefore considered the corresponding auxiliary 
model, in order to obtain this flat sampling. 

Since we want to study the configurations at a given 
energy, the auxiliary model has two terms: the first one 
is quadratic, constraining the energy ( fl9| ) around a given 
Eq-, and the second one is the number of spins not aligned 
with their local field, i.e. the number of spins that can 
flip without thermal activation: 
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E„ 



Pi -J ^ s i s j ~ X! hlSi - E o\ + 



<%,j> 



PzYsSi-SiHi) , 



(20) 



FIG. 18. Energy versus time for the dynamics of the RFIM 
at inverse temperatures 3, 4 and 5. The energy decreases 
slower for larger f3 because the dynamics is activated. For 
one given sample E(t) would be a succession of plateaus, here 
the curves correspond to an average over 64 realizations of 
the random field. 



where is the Heavyside function, and Hi 



<j,i> a 3 



s.j + hi is the local field at site i. The two 



auxiliary inverse temperatures Pi and P% are used to 
perform an annealing, starting from a random initial 
configuration of high energy. We use a single-spin flip 
Metropolis algorithm, accepting the moves with prob- 
ability min(l, exp(— AE aux )). E$ is taken negative and 
slightly larger than the ground state mean energy Eos = 
— 6 JN, in order to look at configurations having the same 
energy than the long time dynamics configurations (the 
evolution of the energy during the dynamical evolution 
is displayed in Fig. p~8| ) . 

The first, simplest observable to look at is the dis- 
tribution of the magnetization, P e ^ w (M), averaged over 
realizations of disorder. Typical results are displayed in 
Fig. |l9], for the P(M) obtained dynamically or with a 
sampling of the blocked configurations, for a system of 
20 3 spins, and for E = —5AJN and -5.2J7V (we have 
also simulated systems of 30 3 spins). While the dynam- 
ical Pdyn(M) consists in a single peak around M = 
(getting narrower for larger system sizes), the distribu- 
tion P e dw{M) is clearly bimodal, presenting two peaks 
around two symmetric values of the magnetization. The 
precise values of the peaks depend on E (decreasing for 
increasing E Q : it is clear that, for too high values of E , 
the configurations are less magnetized), and their width 
depend on system size, getting narrower for larger sys- 
tems. 
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It appears therefore that the configurations dominat- 
ing Edwards' distribution are magnetized. (Note that 
a similar result was obtained in |38j for the low-energy 
metastable states for a ferromagnet on random thin 
graphs). At variance with the previously studied models 
(KA and Tetris), Edwards' distribution is therefore un- 
able to describe the typical configurations obtained dy- 
namically. 

The example of the RFIM clearly underlines the dif- 
ference with the inherent structure construction of Kob 
et al. On the one hand, we have seen that Edwards' 
measure does not reproduce the dynamic results. On the 
other hand, it is clear that different thermal histories, 
just as those considered in [^| . would yield domain con- 
figurations with essentially the same shapes if the end 
energy is the same, and zero magnetization, contrary to 
the typical configurations at that energy. 



0.08 




FIG. 19. Histograms of the magnetization of the visited 
configurations during the low temperature dynamics (full line: 
Pdyn(M) for L = 20, = 4, i.e. the energy per site of the sys- 
tem is between —5.2 and —5.4; the finite width of the peak 
comes from the relatively small size of the sample) and for 
Edwards' measure at Eq — —5.2 and —5.4 (dotted, dashed, 
long-dashed and dot-dashed lines: P e d w (M) for L = 20 and 
02 — 4 and 6; Eq — —5.4 corresponds to the peaks at higher 
M). The dynamics samples configurations with low magne- 
tization, while the typical blocked configurations are magne- 
tized. 
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V. 'CHAOTICITY' PROPERTIES 

We have shown two models for which Edwards' con- 
struction gives a good approximation, and one for which 
it does not. What is the distinguishing feature between 
them? A distinction one can make, suggested by glass 
theory [^9|- ^l| , p0| , is obtained by studying their 'chaotic- 
ity' properties as follows: after aging for a time t w , two 
copies (clones) are made of the system, and allowed to 
evolve subsequently with different realizations of the ran- 
domness in the updating procedure. The question is then 
whether the trajectories diverge or not. Note that for 
this criterion to make sense, it should always be applied 
at non-zero (though weak) tapping or shearing. 

The results summarized below seem to indicate that 
the condition of chaoticity is necessary. It is however not 
sufficient: Bouchaud's 'trap model' 15 is chaotic but its 
fluctuation-dissipation properties are not directly related 
to the density of states Q . 

For the three models discussed in this paper, we have 
measured the normalized average overlap Q tw (t) defined 
as 

• for the KA model 

Qtjt) = d X^Wnf^t) - P (t) • Pclone (t)) 



and in particular for the RFIM, where we have checked 
that the two copies retain a finite overlap at all the times 
studied. 



O 




(21) 



FIG. 20. Mean overlap Qt w (t) between two clones in the 
Tetris model: the two clones are separated at t m and evolve 
subsequently with different noises. Qt w (t) always decreases 
to zero (the slower the larger t w ), showing that the clones 
always diverge. 



where m (resp. nf one ) is 1 if there is a particle 
on the site i for the original model (resp. for the 
clone), and if the site is empty; 



for the Tetris model 



,(*) 



(Pclone(t) ~ B^ll) 



> (22) a 



where Aij{t) is a function that gives one if on the 
site the two copies present the same particles 
and zero otherwise. We have four types of particles 
in the system and this is the reason for the factors 
4 in 



• for the RFIM 

Qtjt) 



clone 



(*)> 



(23) 
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FIG. 21. Same as Fig. M for the KA model. 



VI. CONCLUSIONS 



The brackets indicate an average over different realiza- 
tions of the randomness. 

Fig. and |l] report the results for Qt w (t) for sev- 
eral values of t w , respectively for the Tetris and the KA 
model. It is clear that in both cases the two copies of 
the system always tend to diverge. This is at variance 
with what happens in the domain growth dynamics El , 



In summary, we have proposed a simple and system- 
atic procedure to construct a flat sampling of the 'blocked 
configurations', i.e. to calculate averages with Edwards' 
measure. We have shown, for two representative finite- 
dimensional models, that this measure gives different re- 
sults than the equilibrium measure, and is able to repro- 
duce the dynamical sampling of the out-of-equilibrium 
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compaction dynamics for various observables. The con- 
nection of Edwards' ensemble with the dynamical FDT 
temperature immediately suggests experiments to check 
the validity of these ideas, for example by studying dif- 
fusion and mobility of different tracer particles within 
driven granular media. 

At present, the correspondence between Edwards' dis- 
tribution and the long-time dynamics is at best checked 
but does not follow from any known principle. Now that 
several concrete examples have lent credibility to Ed- 
wards' construction, an effort to understand why it does 
in some cases work and what is its range of validity has 
become worthwhile. 

There remains the question of generalizing Edwards' 
measures in two directions: by considering a fraction 
e > of mobile particles, and by conditioning the flat 
measure to more macroscopic observables in addition to 
density and energy. 
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